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The  influence  of  microcracks  in  the  annular  groove  region  of  a  kinetic  energy  pene- 
trator  is  examined  through  a  combined  theoretical,  experimental,  and  computational 
study.  Microcracks  resulting  from  plastic  deformation  are  quantified  for  a  typical 
penetrator  made  from  tungsten  alloy  (WA).  The  effect  of  existing  and  proposed  an¬ 
nular  groove  geometries  on  the  stress  field  near  the  penetrator  surface  is  determined 
using  finite  element  techniques,  and  the  effect  of  the  various  stress  fields  on  the 
critical  flaw  size  is  examined.  Recommendations  of  microcrack  size  and  groove  ge¬ 
ometry  are  made  to  improve  the  fracture  resistance  of  conventional  WA  long  rod 
penetrators. 


Introduction 

Fracture  of  kinetic  energy  (KE)  penetrators  has  been  observed  in  both  medium- 
and  large-caliber  penetrators.  A  common  feature  between  these  penetrators  is  the  annu¬ 
lar  grooves  used  to  transfer  shear  from  the  sabot  to  the  penetrator  during  launch.  This 
relatively  sharp  feature  acts  as  a  stress  concentration  and  can  lead  to  Mode  I  fracture  in 
the  penetrator  under  certain  loading  conditions  from  either  the  cannon  launch  or  the 
terminal  ballistic  engagement  with  the  annored  vehicle.  Understanding  and  predicting 
this  mode  of  failure  in  high-density  metals  is  critical  to  the  proper  design  of  effective 
long  rod  penetrators.  This  paper  examines  the  state  of  stress  and  strain  at  the  root  of 
several  groove  configurations  for  a  KE  penetrator  made  from  a  tungsten-sintered  alloy. 
The  existence  of  microcracks  near  the  surface  of  the  material,  resulting  from  an  applied 
load,  is  first  examined  through  experiments.  The  conditions  necessary  for  the  micro¬ 
cracks  to  cause  macroscopic  fracture  in  this  class  of  material  is  then  considered.  Fi¬ 
nally,  the  results  from  a  series  of  numerical  simulations  that  model  several  groove  con¬ 
figurations  are  presented  to  identify  the  significance  of  microcracking  on  overall  pene¬ 
trator  structural  integrity. 


Beam  Impact  Experiments 


A  key  parameter  for  the  modeling  presented  in  this  paper  is  the  degree  of  mi¬ 
crocracking,  if  any,  that  exists  in  the  groove  region  of  the  penetrator  material.  To  ad¬ 
dress  this  issue,  a  series  of  experiments  were  performed  to  determine  if  microcracks  do 
in  fact  exist  prior  to  the  formation  of  a  large,  unstable  crack,  and  how  far  beneath  the 
material  surface  the  microcracks  extend.  The  experimental  configuration  used  was  a 
dynamic  beam  impact  test,  depicted  in  Fig  1.  The  beam-shaped  WA  test  specimen  was 
struck  with  a  WA  projectile  of  similar  composition  to  induce  a  dynamic  Mode  I  frac¬ 
ture  event,  initiating  on  the  side  opposite  that  struck  by  the  projectile.  During  dynamic 
beam  flexure  but  prior  to  failure  (i.e.,  unstable  macroscopic  crack  growth),  small  mi¬ 
crocracks  that  either  pre-exist  or  form  during  dynamic  loading  will  open  and  possibly 
remain  open  for  posttest  inspection.  It  is  through  posttest  optical  inspection  of  the  re¬ 
gion  near  the  fracture  surface  that  the  existence  of  microcracks  can  be  verified,  and 
their  depth  determined. 


FIGURE  1 .  Schematic  depiction  of  the  beam  impact  experiments. 

The  experiments  were  performed  using  a  93%  (by  weight)  tungsten,  5.6% 
nickel,  and  1.4%  iron  WHA.  The  bulk  density  was  17.76  g/ctn3.  The  original  material 
was  in  the  form  of  a  bar  that  had  been  swaged  to  20%  reduction  in  area  and  aged.  The 
quasi-static  ultimate  tensile  true  stress  was  1 .425  GPa,  with  a  logarithmic  strain  of 
12.5%,  and  the  quasi-static  yield  stress  was  1.19  GPa  (these  were  furnished  by  the  ma¬ 
terial  supplier).  The  dilatational  and  shear  wave  speeds  were  measured  to  be  5 163  and 
2824  m/s  using  an  ultrasound  technique.  The  material  elastic  modulus  and  Poisson’s 
ratio  were  365  GPa  and  0.29,  respectively.  These  were  determined  from  the  measured 
wave  speeds.  Beam  specimens  were  cut  from  a  large  WHA  bar  using  conventional  ma¬ 
chining  techniques  similar  to  those  used  to  cut  grooves  in  KE  penetrators.  The  beam 
specimen  had  a  length  of  127  mm  and  a  square  cross  section  with  height  and  width  of  9 
mm  each.  The  length  and  diameter  of  the  projectile  were  44.5  and  8.9  mm,  respec¬ 
tively.  The  impact  speed  was  55  m/s. 

The  rear  surface  region  of  the  impacted  beam  specimens  near  the  fracture  plane 
was  examined  with  a  scanning  electron  microscope  to  confirm  the  presence  of  micro- 


cracks  and  to  estimate  microcracking  depth.  A  typical  image  of  this  region  is  shown  in 
Fig  2a.  The  fracture  plane  is  located  just  beyond  the  bottom  of  the  image.  Significant 
microcracking  is  evident,  with  the  microcracks  oriented  perpendicular  to  the  beam  lon¬ 
gitudinal  axis.  The  maximum  crack  opening  was  determined  to  be  approximately 
1 .7  pm.  Images  of  the  beam  surface  taken  far  from  the  region  of  fracture  (not  shown) 
did  not  show  any  evidence  of  microcracking,  suggesting  that  the  large  tensile  stress 
near  the  fracture  plane  of  the  beam  was  responsible  for  either  opening  existing  micro¬ 
cracks  or  causing  them  to  form. 

The  depth  of  the  microcracking  was  estimated  from  additional  micrographs 
taken  along  the  region  where  the  macroscopic  fracture  event  initiated.  Fig  2b  shows  a 
typical  surface  topography  of  this  region.  The  beam  specimen  was  rotated  into  the 
plane  of  the  image  by  40°  to  obtain  the  necessary  perspective  to  estimate  microcracking 
depth.  The  circled  area  of  Fig  2b  shows  the  remnants  of  a  microcrack.  As  can  be  seen, 
the  crevasse  has  a  depth  of  one  tungsten  grain,  or  -40  pm.  The  formation  of  micro¬ 
cracks  in  a  tungsten-sintered  alloy  during  tensile  loading  was  discussed  in  detail  by 
Weerasooriya  [1],  who  found  the  grain-to-grain  contact  to  be  the  weakest  portion  of  the 
alloy,  and  therefore  the  first  to  fail.  An  identical  mechanism  is  likely  responsible  for 
the  microcracking  observed  in  the  recovered  beam  impact  specimens. 


(a)  (b) 

FIGURE  2.  Scanning  electron  microscope  image  of  the  beam  specimen,  (a)  Surface  near  the  fracture 
plane.  Actual  fracture  plane  is  located  just  beyond  the  bottom  of  the  image,  (b)  Oblique  view 
of  a  microcrack  interacting  with  the  main  fracture  plane.  Circled  region  identifies  the  microcrack 

of  interest. 


The  beam  impact  experiments  verily  that  microcracks  can  exist  in  tungsten- 
sintered  alloys  when  subjected  to  tensile  loads  below  the  catastrophic  limit.  Further¬ 
more,  microcracking  appears  to  penetrate  to  a  depth  of  one  grain  diameter  below  the 


machined  surface.  This  provides  the  estimates  of  microcrack  length  needed  for  the 
elastic -plastic  fracture  mechanics  and  finite  element  analyses  which  are  presented  next. 


Elastic-Plastic  Fracture  Mechanics  Considerations 

The  observations  from  the  beam  impact  experiments  confirmed  the  presence  of 
microcracks,  having  a  depth  of  approximately  one  tungsten  grain  diameter  (40  //in)  be¬ 
low  the  surface.  In  the  context  of  a  KE  penetrator,  these  microcracks  will  exist  in  the 
groove  region  of  the  penetrator,  which  is  subject  to  macroscopic  elastic-plastic  defor¬ 
mation.  Hutchinson  [2]  and  Rice  and  Rosengren  [3]  each  developed  crack  tip  solutions 
for  materials  capable  of  plastic  behavior.  Introducing  plasticity  eliminated  the  possibil¬ 
ity  for  using  linear  elastic  fracture  mechanics  (LEFM)  in  the  direct  vicinity  of  the  crack 
tip.  Far  from  the  crack  tip  (but  still  in  the  crack  tip  region)  LEFM  does  apply.  How¬ 
ever,  the  region  near  the  crack  tip  is  characterized  by  the  plastic  solutions  of  Hutchin¬ 
son,  Rice,  and  Rosengren,  and  is  referred  to  as  the  HRR  field.  The  stress  field  ahead  of 
a  sharp  crack  in  an  elastic-plastic  material  is  shown  in  Fig  3.  The  decrease  in  stress 
very  close  to  the  crack  tip,  denoted  as  RJ,  is  due  to  crack  tip  blunting  [4], 


FIGURE  3.  Schematic  depiction  of  the  opening  stress  field  ahead  of  a  sharp  crack. 

Application  of  LEFM  and  the  HRR  field  to  determine  the  fracture  initiation  cri¬ 
terion  for  low  temperature  mild  steel  was  perfonned  by  Ritchie,  Knott,  and  Rice  [5], 

At  low  temperatures,  mild  steel  (BCC  crystal  structure)  fractures  by  brittle  grain  cleav¬ 
age  [6],  This  crystal  structure  and  mode  of  failure  are  similar  to  those  of  the  tungsten- 
sintered  alloy  examined  in  this  study.  Ritchie,  Knott  and  Rice  found  that  a  critical 
stress  must  be  exceeded  at  a  certain  distance  ahead  of  the  crack  tip  for  fracture  initiation 
to  occur.  The  critical  stress  and  distance  are  denoted  in  Fig  3  as  rc  and  ac ,  respect- 


tively.  For  brittle  fracture  of  mild  steel,  rc  was  found  to  be  two  grain  diameters.  It  was 
hypothesized  by  Ritchie,  Knott,  and  Rice  that  the  driving  force  necessary  for  grain 
cleavage  ( a c )  must  extend  beyond  the  first  grain  due  to  the  random  grain  orientation  of 

the  material  and  the  preferential  cleavage  directions  of  BCC  grains.  Should  the  first 
grain  cleave,  the  crack  tip  would  likely  be  in  front  of  a  grain,  oriented  in  a  way  that 
would  be  less  favorable  to  cleavage.  Hence,  an  overabundance  of  driving  energy  is 
needed  for  macroscopic  fracture  to  occur.  A  similar  argument  for  unstable  crack  initia¬ 
tion  is  proposed  for  the  tungsten-sintered  alloy  examined  here.  In  what  follows,  a  finite 
element  technique  is  used  to  examine  both  the  crack  tip  driving  energy  for  an  elastic- 
plastic  material  (the  J-integral)  as  well  as  the  extent  to  which  the  stress  field  penetrates 
below  the  material  surface  for  various  groove  configurations  and  microcrack  depths. 
The  intent  here  is  not  to  quantify  the  fracture  criterion,  but  rather  to  examine  the  rela¬ 
tive  effects  of  various  groove  configurations  and  microcrack  depths. 


Finite  Element  Analyses 

The  state  of  stress  and  strain  in  the  groove  region  of  a  KE  penetrator  was  exam¬ 
ined  using  the  finite  element  code  ANSYS.  For  large-caliber  penetrators,  the  grooves 
typically  have  a  7°  driving  land  and  a  45°  clearance  flank  with  blend  radii  to  the  minor 
diameter,  as  shown  in  Fig  4.  Because  these  radii  act  as  stress  risers  under  bending 
loads,  changing  to  a  triradial  blend  between  the  7  and  45°  surfaces  has  been  demon¬ 
strated  to  decrease  bending  stresses  by  20-25%;  thus,  possibly  decreasing  the  energy 
required  to  initiate  an  existing  flaw  into  a  propagating  crack.  An  initial,  baseline  nu¬ 
merical  simulation  was  performed  using  grooves  with  a  single  and  triradial  blend  on  a 
penetrator  subjected  to  quasi-static  three-point  bend  loading.  To  simplify  the  finite 
element  model,  a  harmonic,  axisymmetric  element  was  used.  While  this  reduced  prob¬ 
lem  size  and  permitted  adequate  mesh  resolution,  the  harmonic  element  limited  the 
analysis  to  linear  elastic  response.  This  initial  analysis  indicated  that  the  triradial  pro¬ 
file  reduced  maximum  bending  stresses  by  17%,  but  the  stress  profiles  were  the  same 
for  locations  100  pm  or  greater  below  the  surface.  Furthennore,  the  stresses  were  sev¬ 
eral  times  the  yield  stress,  indicating  the  need  to  introduce  material  plasticity  into  the 
simulations. 

Subsequent  simulations  utilized  modeling  techniques  that  permitted  material 
plasticity  and  the  inclusion  of  sharp  cracks.  A  plane  stress  approach  was  used  and  the 
geometry  was  mirrored  about  the  axis  of  symmetry.  By  setting  element  thickness  to 
maintain  cross-sectional  inertia  of  the  “circular”  penetrator,  deflections  resulting  from 
bending  loads  would  be  accurate  in  the  “rectangular”  penetrator  model.  Finally,  be¬ 
cause  the  distance  to  the  neutral  axis  was  also  preserved,  stresses  would  also  be  main¬ 
tained.  A  detailed  analysis  of  this  modeling  technique  is  presented  in  Sorensen  [7]. 


Modified 


FIGURE  4.  Groove  profile  illustrating  the  7  and  45°  flanks  and  a  superposition  of  the  original  and 
modified  groove  profiles  illustrating  crack  locations. 

Before  jumping  straight  into  the  complexity  of  elastic-plastic  numerical  simula¬ 
tions  with  pre-existing  cracks,  a  set  of  linear  elastic  simulations  with  cracks  of  a  given 
size  were  first  perfonned.  This  provided  an  understanding  of  the  structural  response  to 
specific  flaws,  or  more  specifically,  to  flaw  location  and  size.  For  the  elastic  analysis, 
infinitely  sharp  cracks  were  introduced  into  three  locations  of  each  groove  profile,  as 
shown  in  Fig  4.  These  locations  were  selected  because  they  coincide  with  the  maxi¬ 
mum  stress  locations  in  both  groove  profiles:  A  and  B  for  the  original  groove  and  C  for 
the  modified  groove.  In  each  case,  the  crack  was  normal  to  the  surface  and  at  the  same 
axial  locations  for  either  groove  profile. 

Results  from  these  simulations  are  shown  in  Fig  5.  In  Fig  5a,  the  J-integral  is 
plotted  as  a  function  of  crack  length  and  illustrates  global  response  to  crack  length.  In 
Fig  5b,  the  J-integrals  are  nonnalized  by  the  J-integral  from  crack  C  for  the  modified 
profile  and  show  relative  response  to  crack  length  and  location.  Note  that  J  is  propor- 
tionnal  to  K]  ,  where  Kj  is  the  Mode  I  stress  intensity  factor  because  the  analysis  is  lim¬ 
ited  to  linear  elastic  material  behavior  [8].  The  primary  observation  is  that,  as  expected, 
the  J-integral  increases  with  crack  length.  Secondly,  depending  on  location  and  extent 
of  stress  risers,  the  relative  magnitude  of  maximum  stress  intensity  varies  by  up  to  30% 
between  the  two  grooves.  However,  if  critical  flaw  size  is  in  excess  of  40  pm  (i.e.,  one 
tungsten  grain  diameter),  this  difference  decreases  to  15%.  Therefore,  from  a  LEFM 
analysis,  the  difference  between  the  two  profiles  may  be  irrelevant,  particularly  if  there 
is  an  overabundance  of  driving  energy. 

To  begin  to  understand  the  effects  of  plasticity,  two  limited  numerical  analyses 
were  performed  using  a  bilinear  elastic-plastic  model.  The  first  analysis  examined  the 
plastic  response  of  the  penetrator  without  the  imposed  flaw.  The  results  showed  that 
for  the  applied  loading,  maximum  tensile  plastic  strains  of  4.5%  and  3.0%  were  present 
for  the  original  and  modified  grooves,  respectively.  However,  while  the  plastic  strains 
in  the  original  groove  were  50%  higher,  the  extent  of  the  plastic  strain  field  in  excess  of 
3.0%  (i.e.,  the  size  of  the  plastic  zone)  was  smaller  for  the  modified  groove. 


Crack  Length  (pm)  Crack  Length  (pm) 


(b) 

FIGURE  5.  Results  from  linear-elastic  fracture  analysis. 

The  second  numerical  analysis  examined  the  effects  of  crack  size  with 
plasticity  present.  Here,  only  the  crack  location  in  the  maximum  stress  riser  was  exam¬ 
ined  for  crack  sizes  up  to  100  pm  (locations  A  and  C  for  the  original  and  modified 
grooves,  respectively).  Furthermore,  the  applied  load  had  to  be  reduced  50%  to  avoid 
instabilities  during  solution.  Table  1  presents  results  in  terms  of  the  maximum  effective 
plastic  strain  at  the  crack  tip  and  the  extent  of  the  plastic  strain  field  beyond  the  crack 
tip  as  a  function  of  crack  length.  The  results  in  Table  1  indicate  that  there  is  a  demarca¬ 
tion  for  plastic  strain  and  plastic  zone  extent  with  respect  to  crack  length  between  30  to 
50  pm.  For  crack  lengths  greater  than  50pm,  plastic  strain  and  its  extent  are  moving 
toward  an  asymptote,  and  the  differences  between  the  two  grooves  are  decreasing.  Fur¬ 
thermore,  the  extent  of  the  plastic  zone  is  in  excess  of  the  average  grain  size.  The  re¬ 
sults  do  show,  however,  that  for  small  cracks,  the  maximum  strain  and  plastic  zone  size 
are  slightly  smaller  for  the  modified  groove  design.  Hence,  if  the  fracture  initiation  cri¬ 
terion  of  Ritchie,  Knott,  and  Rice  discussed  earlier  is  true  for  tungsten-sintered  materi¬ 
als,  use  of  the  modified  groove  geometry  may  result  in  a  penetrator  design  more  resis¬ 
tant  to  fracture.  However  the  differences  for  the  larger  crack  sizes  may  be  sufficiently 
small  to  have  negligible  influence  on  fracture  initiation. 


TABLE  1.  Plastic  Strain  Field  Parameters  for  Various  Crack  Sizes 


Crack 

Length 

(pm) 

Max.  Effective  Plastic  Strain 
(%) 

Plastic  Strain  Field  Size 
(pm) 

Original 

Groove 

Modified 

Groove 

Original 

Groove 

Modified 

Groove 

25 

69 

60 

35 

25 

50 

92 

87 

75 

55 

75 

105 

106 

100 

90 

100 

114 

121 

110 

115 

Conclusions 

The  existence  of  microcracks  in  the  high-stress  region  surrounding  a  macro¬ 
scopic  fracture  in  a  93%  tungsten- sintered  alloy  was  verified  through  an  experimental 
technique.  Their  presence  suggests  that  a  fracture  initiation  criterion  that  includes  not 
only  maximum  stress  or  strain  near  the  crack  tip,  but  also  the  size  of  the  region  encom¬ 
passed  by  high  stresses  or  strains  may  apply.  Elastic-plastic  numerical  simulations  of 
various  groove  configurations  with  pre-existing  microcracks  indicated  a  modified 
groove  geometry  would  slightly  decrease  the  maximum  strain  and  size  of  the  plastic 
zone;  however,  the  decrease  may  not  be  sufficient  to  affect  fracture  initiation. 
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